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Using particle simulations of the uniform shear flow of a rough dilute granular gas, we 
show that the translational and rotational velocities are strongly correlated in direction, but 
there is no orientational correlation-induced singularity at perfectly smooth (f3 = —1) and 
rough (f3 = 1) limits for elastic collisions (e = 1); both the translational and rotational ve- 
locity distribution functions remain close to a Gaussian for these two limiting cases. Away 
from these two limits, the orientational as well as spatial velocity correlations are respon- 
sible for the emergence of non-Gaussian high velocity tails. The tails of both distribution 
functions follow stretched exponentials, with the exponents depending on normal (e) and 
tangential (/3) restitution coefficients. 

PACS numbers: 45.70.Mg, 47.45.-n, 45.70.-n 

Under external forcing (e.g. shearing/vibration), the granular materials, a collection of macro- 
scopic solid particles, can flow like a gas or a liquid. The rapid flow of granular gases (dilute 
limit) has been extensively studied using kinetic theory that takes into account the dissi- 

pative nature of particle interactions. Most theories of granular gases assume that the particles are 
smooth which is a simplification of real particles which are always rough, giving rise to surface 
friction, and the rotational motion is important to deal with such rough, frictional particles. Prior 
literature M, \% a [60 suggests that the rotational motion should not be neglected for a realistic 
modeling of the dynamics and pattern formation in a granular gas even in the Boltzmann (dilute) 
limit. For example, the molecular dynamics (MD) simulations [6] have elucidated the crucial role 
of friction on pattern formation in oscillated granular layers via a comparison with experimental 
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results. Despite the importance of frictional interactions, only a small-body of work exists on the 
modeling of rough granular gas [|2|,|3|,|4j,|50. 

To develop constitutive models of rough granular gases, a systematic study of correlations 
and the distribution functions of both 'translational' and 'rotational' velocities is of fundamental 
interest. While the deviation of translational velocity distribution functions (VDF) from a Gaussian 
has been extensively studied (in terms of stretched exponential or power-law tails) using theory ||7|], 
simulation |8|] and experiment [9], similar results on 'rotational' VDFs are very scarce. 

For a 'rough' granular gas, one needs to probe possible 'orientational'/'directional' correlations 
between translation and rotation, in addition to standard density and velocity correlations. It has 
been recently found |5Q that such orientational correlations are strong and the limit of smooth 
granular gas is singular in a freely cooling granular gas. The last result readily raises doubts 
about the validity of the hydrodynamic theories Q4£| that are obtained via perturbative expansions 
around the smooth-particle limit. Also, it is of interest to ascertain the impact of such orientational 
correlation on VDFs, especially in two limits of perfectly smooth and rough particles around 
which the perturbative expansions are sought. The above issues are investigated in this paper using 
MD simulations of the uniform shear flow (which is a prototype non-equilibrium steady-state to 
develop constitutive relations [4]) of a dilute rough granular gas. 

We consider a mono-disperse system of rough, inelastic spheres of diameter d, mass m, and the 
moment of inertia X, interacting via purely repulsive potential. Let us denote the pre-collisional 
translational and rotational velocities of particle i by and uji, respectively, and the corresponding 
post-collisional velocities are denoted by the primed symbols, c[ and uj[. The pre-collisional 
relative velocity at contact, gij, between particle i and j is g^ = c„ — ((i/2)kx(w, + u>j), where 

= Cj — Cj is the relative translational velocity between particle i and j, and k is the unit vector 
directed from the center of particle j to that of particle i. Neglecting Coulomb friction in the dilute 
limit, we have the following collision model[|2j]: 

k-g-j = -e(k.gy) and kxg^ = -/?(kxg^), 

characterized by two parameters: the normal restitution coefficient, < e < 1, and the tan- 
gential restitution coefficient, — 1 < @ < 1. The former is an indicator of the inelasticity of 
particle-collisions and the latter characterizes its surface roughness. For collisions between per- 
fectly smooth particles @ — — 1, with increasing value of (3 being an indicator of the increasing 
degrees of particle surface friction. The value of (3 = represent the case for which the particle 
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surface friction and inelasticity are sufficient to eliminate the post-collisional tangential relative 
velocities. For < (5 < 1, the spin-reversal occurs after collision[2], and the case of /3 = 1 
corresponds to collisions between perfectly rough particles. 

We have used an event-driven algorithm [10] to simulate the uniform shear flow, characterized 
by a linear velocity profile. The simulation box is a cube, with Lees-Edwards boundary condi- 
tion nilll across two moving boundaries along y-direction, and the periodic boundary conditions 
along x- and z-directions, with x being the direction of flow. The positions of particles are initial- 
ized in the simulation box, with their translational and rotational velocities being chosen randomly 
from a Gaussian distribution. The data are accumulated once the system has reached a statistical 
steady-state condition which is monitored from the temporal evolution of the system's kinetic en- 
ergy. The coarse-grained 'translational' fluctuation kinetic energy (i.e. the standard granular tem- 
perature ), T, and the 'rotational' fluctuation kinetic energy, 9, are defined as: T(x, t) = (C • C)/3 
and 6>(x, t) = (I/3m)(Cl-Cl), respectively, where the angular bracket denotes a suitable averaging. 
Here C = c — (c) is the translational 'peculiar' velocity which measures the deviation of the in- 
stantaneous particle velocity from the local mean velocity ((c)), and = oj — (u>) is its rotational 
counterpart. All results are presented for a very dilute system (Boltzmann limit), with a volume 
fraction of = 0.01. (The effect of density, as well as Coulomb friction, is a non-trivial issue 
which will be taken up in a separate publication.) The number of particles was fixed to N = 8000, 
and the robustness of results was checked by using iV = 4000 and 16000. 
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FIG. 1: (color online) (a) Variation of A(t) = cos 2 \P with time for different values of e with (3 = 0, 
4> = 0.01 and N = 8000. Left and right insets show the distributions of cos 2 4* and cos 4 1 , respectively, (b) 
Variation of (cos 2 \&) with for different e. Larger symbols (triangle and hatched-circle) at j3 = for each 
e correspond to simulations with N = 4000 and 16000, respectively. 
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It is known that the translational and rotational fluctuating velocities are uncorrelated in a 
molecular gas, but they have been shown to be correlated (in direction) in a freely cooling granular 
gas B50. This orientational/directional correlation between translational and rotational velocities is 
quantified in terms of the mean square of the cosine of the angle, between C = c — (c) and 

= uj — {uj)\ 

i=i x 1 ii i=i 

In fig. dja) we have plotted the temporal variation of A(£) (main panel) for two values of normal 
restitution coefficients (e = 1, 0.5), with the tangential restitution coefficient being set to (3 — 0; 
the corresponding probability distribution of cos 2 \I/ is shown in the left inset. (The probability 
distribution of cos P(cos is symmetric about its zero mean for all e, but its width becomes 
narrower with decreasing e, see the right inset.) From the main panel and the left inset, we find 
that even for e = 1 the mean value of A is different from 1/3 (for a molecular gas), signaling 
the presence of orientational/directional correlation; decreasing the value of e to 0.5 decreases 
its value to (A) ~ 0.26, thus enhancing orientational correlation significantly. The variation of 
the temporal- average of A(£) with particle roughness, (3, is shown in fig. \Hb); the dot-dash line 
represents the limiting value of 1/3 for a molecular gas. Note that the data points for e = 1 
(thick blue dashed line) and e = 0.99 almost overlap with each other. For any e, the orientational 
correlation is maximum at (3 ~ and it decreases monotonically as we approach the perfectly 
smooth (/3 = —1) and perfectly rough {(3 = 1) limits. This latter observation is in contrast to 
that in a freely cooling dilute granular gas [|50 for which (A) varies non-monotonically with (3 for 
— 1 < (3 < and < (3 < 1 . Another difference with freely cooling gas is that the magnitude of 
(A) is much larger in shear flow. It must be noted that even though the translation and rotation are 
decoupled at (3 = — 1 (independent of the value of e), the smooth limit is singular for any e ^ 1 in 
shear flow. However, there is no orientational correlation-induced singularity at both the perfectly 
smooth {[3 = —1) and rough {[3 = 1) limits for the limiting case of e = 1. 

Before presenting results on VDFs, we probe mean-field quantities to ascertain the presence 
of any inhomogeneity in our system. Figure [2] shows the probability distributions of the mean 
density (main panel) and the rotational temperature (right upper inset) for = 0.01 and (3 = 0; the 
distribution of mean translational temperature (T) looks similar to that of 9 (not shown). The lower 
right inset shows that the translational and rotational temperatures are un-equally partitioned over 
the whole range of (3 (except at j3 = 1 with e = 1), and the calculated temperature ratio (symbols) 
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FIG. 2: (color online) Probability distributions of mean density (main panel) and rotational temperature 
(right upper inset) for different values of e. Right lower inset shows the variation of temperature ratio, 
6/(9 +T), with p. Two left insets show projected particle snapshots in the xy- and yz-planes at steady 
states for e = 0.5. 

agrees well with theoretical predictions (solid line). The data for mean distributions in Fig. 2 
have been obtained by dividing the simulation box into a number of equal-sized cells (10 3 ) such 
that on average about ten particles occupy each cell and then calculating the instantaneous value 
of any mean field quantity (0, T, 9, (c)) in each cell. It is interesting that even though the mean 
density varies between 0.003 and 0.017 in about 90% cells, the density-distribution remains almost 
identical with decreasing e from 0.99 to 0.5. The projected snapshots of all particles (at steady 
state after 60000 collisions per particle) in the xy- and ?/z-planes, as displayed in two left insets of 
Fig. 2, further suggest that the particles are homogeneously distributed and there is no discernible 
dissipation-induced clustering even at e = 0.5 in our system. 

Now we turn to velocity distribution functions which are known to be strongly affected by the 
presence of correlations [0, M, ISO- To calculate VDFs, we used cell-wise averaging as discussed 
above. Figure 3(a) shows the probability distribution functions of the translational (main panel) 
and rotational (lower inset) fluctuating velocities for (3 = —0.999, while Fig. 3(6) shows the trans- 
lational (lower inset) and rotational (main panel) VDFs for [3 = 0. In each plot, two data sets for 
e = 0.99 and 0.5 have been superimposed, and the black dashed line represents a Gaussian. The 
horizontal axis of each plot has been scaled by a, the standard deviation of the corresponding dis- 
tribution, and the vertical axis scaled such that P(0) = 1. For the perfectly smooth limit {(3 ~ —1) 
in Fig. 3(a), the tails of the translational VDFs deviate from a Gaussian with increasing dissipa- 
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FIG. 3: (color online) (a) Translational (main panel) and rotational (lower inset) VDFs for j3 = —0.999. (b) 
Rotational (main panel) and translational (lower inset) VDFs for (3 = 0. The upper right inset in each panel 
shows the deviation of the low velocity regions from a Gaussian. 

tion; the rotational VDFs remain a Gaussian even for e = 0.5. We did not find any discernible 
difference among the VDFs for all three components of each velocity (Cj and fij, with i = x,y, z) 
and that they follow the same distribution for a given e. For the other extreme of perfectly rough 
limit (P ~ 1, not shown), the translational VDFs follow a similar behavior as that for (3 ~ — 1, 
but the rotational VDFs become non-Gaussian with increasing dissipation. At intermediate values 
of roughness {(3 = 0) as in Fig. 3(6), both the translational and rotational velocity distributions 
deviate from a Gaussian even in the quasi-elastic limit (e ~ 1). It is seen that both P(C X ) and 
P(£l z ) are under-populated for low-velocities and over-populated for high-velocity tails. The up- 
per right insets in Figs. 3(a) and 3(6) show the deviation of low velocities from a Gaussian (i.e. 
AP(x) = P(x) — exp(— x 2 /2), with x = C/a,Q/a). It is noteworthy that the low-velocity re- 
gion of P(Q Z ) is slightly asymmetric for e = 0.5; such asymmetry is absent for £l x and £l y and 
at higher values of e > 0.7. Such asymmetry implies the onset of preferential transport of rota- 
tional velocity fluctuation along the negative ^-direction (i.e. the mean vorticity direction of the 
steady shear flow), but the reason for the emergence of this asymmetry (at large dissipation levels) 
remains unclear with periodic boundary conditions along z. 

Our data on high-velocity tails of P(C) and P(Q) have been fitted with stretched exponential 
functions [7] of the form P(x) ~ exp(— ^ x x ax ), where a x and ^ x (with x = C, Q) are the exponent 
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FIG. 4: (color online) Variation of the exponent for the stretched exponential, a^, with (3; translational (main 
panel, ac) and rotational (inset, an) velocity. Typical error-bars are shown on few data points. 

and pre-factor of the corresponding distribution. Figure 4 shows the variations of ac (main panel) 
and an (inset) with roughness (3 for different e. For e = 1 the tails of both translational and 
rotational VDFs deviate from a Gaussian (except at (3 — ±1 for which ac = 2 = an), and 
the functional-forms of a c and an are asymmetric and symmetric (around /3 = 0), respectively. 
With increasing inelasticity, both ac and a n decrease sharply, and «sj also becomes asymmetric 
around (3 = 0. A least-square fit to our data suggests that ai follows a power-law relation with 
inelasticity: an = 2 - A^l - e 2 ) Bl with i = C,Q, and (A C ,B C ) « (5/8,2/3) at /3 = -1, and 
(A c , B c ) « (3/4, 1) and (Aq, ifo) « (3/8, 7/8) at /5 = 1. Due to the asymmetry of a c and a n 
around (3 = 0, we could not find an universal scaling of ac and ckq with (1 — /5 2 ) at any e. 

Lastly, the variations of spatial translational velocity correlation function, GQ C {r/d) = 
{C X (R)C X (R + r)), for different e are shown in fig. 5 for (3 = (main panel) and /? = —0.999 
(inset)- other velocity components follow a similar behavior. It is seen that Gqc (and related 
correlation length) increases with increasing inelasticity for a given (3, and the effect of rotational 
dissipation {(3) is more prominent in the quasi-elastic limit (e ~ 1). (We have checked that the 
pair correlation function, g(r/d), for any (3 (not shown) is featureless, except that its contact value 
increases slightly with decreasing e.) A general finding is that in the quasi-elastic limit the den- 
sity and spatial velocity correlations are negligible near (3 — ±1, however, the spatial correlations 
for translational velocities emerge with increasing rotational dissipation and becomes strong near 
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FIG. 5: (color online) Spatial velocity correlation: (3 = (main panel) and (3 = —0.999 (inset). 

(3 = 0. These velocity correlations, together with orientational correlations, are responsible for 
non-Gaussian VDFs in a rough granular gas. 

In conclusion, we showed that the translational and rotational velocities are directionally corre- 
lated in a dilute sheared granular gas, but there is no orientational correlation-induced singularity 
at both the perfectly smooth {(3 = — 1) and rough {(3 = 1) limits for e = 1. Even though the 
translational and rotational degrees of freedom are uncorrected at (3 = —1 for any value of e, 
the smooth-limit ((3 — > —1) is singular (i.e. < A >^ 0) for any e ^ 1. Both 'translational' and 
'rotational' VDFs remain close to a Gaussian for these two limiting cases; away from (3 = ±1, 
the orientational correlations and spatial velocity correlations are responsible for the emergence 
of non-Gaussian VDFs. The tails of both VDFs follow stretched exponentials, with the exponents 
depending on two restitution coefficients (e and (3). One immediate consequence of our results 
is the resolution of a doubt [5] that in a dilute sheared granular gas the perturbative expansions 
around the smooth (and also rough) limit is appropriate with Gaussian being the leading-order 
VDF 040. 
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